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Hysteretic superconducting resistive transition in Ba .o7Ko.93Fe 2 As 2 
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We have observed hysteresis in superconducting resistive transition curves of Bao.o7Ko.93Fe2As2 
(T c ~8 K) below about 1 K for in-plane fields. The hysteresis is not observed as the field is tilted 
away from the ab plane by 20° or more. The temperature and angle dependences of the upper critical 
field indicate a strong paramagnetic effect for in-plane fields. We suggest that the hysteresis can 
be attributed to a first-order superconducting transition due to the paramagnetic effect. Magnetic 
torque data are also shown. 

PACS numbers: 74.70.Xa, 74.25.Dw, 74.25.Op 



I. INTRODUCTION 

Magnetic fields destroy spin-singlet superconductivity 
via two different mechanisms: orbital and spin param- 
agnetic effects^ The former leads to the formation of a 
mixed state and, in the absence of the latter, a second- 
order transition to the normal state at the orbital crit- 
ical field B* 2 (0) = $o/27r£ 2 , where <J> and £ are the 
flux quantum and superconducting coherence length, re- 
spectively. The latter lowers the normal state energy 
due to spin polarization and, in the absence of the for- 
mer, may cause a first-order transition to the normal 
state at the paramagnetic critical field B po = A/ \/2/i b , 
where A and fiB are the superconducting energy gap 
and Bohr magneton, respectively. The Maki parameter 
a = V2-B* 2 (0) / B po describes the relative importance of 
the two effects. The three parameters may be estimated 
from experimental data using the weak-coupling BCS re- 
lations: #* 2 (0) = 0.693|^ 2 |T C , B po (in Tesla) =1.84T C , 
and a = 0.528|B£ 2 |, where B' c2 = dB c2 /dT\ Tc . 

It was predicted that when a > 1 the transition from 
the superconducting to normal state becomes first order 
at low temperatures^ It was subsequently proposed that 
in such cases a novel superconducting state, now called 
the FFLO state, in which an order parameter oscillates in 
real space due to a finite center-of-mass momentum q =/= 
of Cooper paris, may occur between the BCS and normal 
states* 3 -^ Although those predictions were already made 
in mid 1960's, strong experimental evidence for them has 
appeared only recently, for two types of compounds: a 
heavy-fermion superconductor CeCoIng^r— and organic 
superconductors J£r— 

Iron-pnictide superconductors, first discovered by 
Kamihara et al.^ are also good candidates for the obser- 
vation of the first-order transition or the FFLO stated 
They have large upper critical fields B c2 , and clear para- 
magnetic limiting of B c2 (T) curves has indeed been re- 
ported for some of them.— _ — Importantly they are multi- 
band superconductors^ In multi-band superconductors, 
high-field behavior due to the spin paramagnetic effect 



would contain much richer physics than that in single- 
band superconductors, as illustrated by the following 
questions.— The Maki parameter may be estimated for 
each band from band parameters. If a > 1 for only one 
band, does the transition become first order ? The FFLO 
q vector is given by q = gfisH/hvF in the simplest case, 
and hence q varies from band to band. How is a compro- 
mise reached in reality ? 

We report here electrical resistance and magnetic 
torque measurements on Baoo.7Ko.93Fe 2 As 2 . It is a 
multi-band superconductor, where the Fermi surface con- 
sists of three large hole cylinders at the zone center and 
small hole cylinders near the corners i^r— The size of the 
superconducting energy gap varies considerably from FS 
cylinder to cylinderj 29 i 30 We observe hysteresis in resis- 
tive transition curves at low temperatures for in-plane 
fields and suggest that it is due to a first-order supercon- 
ducting transition. 



II. EXPERIMENTS 

High-quality Bao.o7Ko.93Fe 2 As 2 single crystals were 
prepared by a self-flux method^ The Ba-to-K ratio was 
determined from energy-dispersive X-ray analysis (aver- 
age of measurements on three pieces at 10 points per 
each piece). The variation in the Ba content was within 
~ ±0.02 from sample to sample and also from point to 
point. Standard four-contact electrical resistance R mea- 
surements were performed in a dilution refrigerator and 
superconducting magnet. A low- frequency ac current 
(typically / = 13 Hz and I = 100 [iA, corresponding to 
the current density in the order of 1 A/cm 2 ) was applied 
in the ab plane and perpendicular to the applied field 
B. The field angle is measured from the c axis: 8 = 
90° for B || ab. Among five measured samples from the 
same growth batch, the resistive hysteresis was clearly 
observed in two, only just observed in one, and not ob- 
served in the last two. The clearly hysteretic samples, 
called HI and H2 hereafter, have the residual resistivity 
ratio at T = 12 K and Tc of (HI) 62 and 8.2 K and of 
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(H2) 74 and 7.4 K. The slight difference in T c is probably 
due to difference in the composition. The two samples 
showing no hysteresis have similar resistivity ratios: 65 
and 75, respectively. However, the transition widths of 
R{B) curves are very different between the two groups. 
The 10% to 90% widths for B || ab at T < 0.05 K are 
0.95 and 0.64 T for HI and H2, respectively, while those 
are 1.8 and 2.8 T for the samples showing no hysteresis. 
As shown below (Figs. 1 and 2), transition curves of the 
hysteretic samples, especially H2, consist of small steps. 
This probably indicates that the samples are composed 
of a small number of domains which are homogeneous 
inside but have slightly different compositions from each 
other and that the steps correspond to respective B^'s of 
the domains. On the other hand, transition curves of the 
samples showing no hysteresis are much broader and fea- 
tureless. It seems that the intrinsic hysteresis is blurred 
by larger compositional inhomogeneity. Magnetic torque 
measurements were performed on small pieces cut from 
samples HI and H2 using piezoresistive microcantilevers. 



III. RESULTS AND DISCUSSION 

Figure 1 (a) shows superconducting resistive transitions 
in sample HI at selected temperatures (except for the top 
curve). At T — 0.03 K, the hysteresis is visible between 
B ~11.8 T, where a finite resistance appears, and ~13.5 
T. There is a kink at ^12.5 T, above which the resistance 
rises steeply. The zero-field transition curve (top curve) 
also shows a similar kink and low-temperature tail below 
it, which might suggest the existence of low-T c and hence 
low-i? C 2 regions in the sample. Five times difference in 
the field sweep rate (0.5 vs 0.1 T/min) leads to no es- 
sential difference in the R(B) curves (compare the two 
lowest curves). Four times difference in the measuring 
current (0.05 vs 0.2 mA) leads to no essential difference 
except for a slight non-Ohmic dependence (inset). These 
observations seem to exclude possible extrinsic origins of 
the hysteresis such as history-dependent inhomogeneous 
current path or induced current during field ramping, 
which would largely be affected by the current density or 
sweep rate. As the temperature is raised, the hysteresis 
is still visible at 0.83 K but not at 1.2 K. Figure 1(b) 
shows R(B) curves for different field orientations. As the 
field is tilted from the ab plane, the hysteresis is visible 
at 9 = 80° but not at 70°. 

Resistance measurements on sample H2 give consis- 
tent results. The resistive hysteresis is observed approx- 
imately up to T = 1 K (Fig. 2) and for field orientations 
within ~ 10° of the ab plane (data not shown). Although 
the transition width is narrower than that in sample HI, 
transition consists of several steps as indicated by clear 
peaks in the derivative dR/dB curves (Fig. 2, bottom 
curve). As noted above, this can be interpreted as an in- 
dication that the sample consists of domains with slightly 
different compositions. 

Figure 3(a) shows the temperature dependences of B C 2 
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FIG. 1. (color online), (a) The top curve shows a zero-field 
superconducting transition in sample HI as a function of tem- 
perature (top axis). The others show resistive transitions in 
fields parallel to the ab plane at selected temperatures as 
a function of field (bottom axis). The curves are vertically 
shifted for clarity. Hysteresis is visible at T — 0.03 and 0.83 
K but not at 1.2 K. The two lowest curves were obtained with 
five times different field sweep rates. A 50% criterion for the 
determination of B C 2 is explained for the upper T — 0.03 K 
curve. The inset shows the current dependence of the transi- 
tion curves, (b) Resistive transition curves for selected field 
angles. The hysteresis is visible at 6 — 90, 87, and 80° but 
not at 70° 



in sample HI. Linear fitting to three data points closest 
to T c gives the initial slopes of B' c2 — -5.3 and -1.1 T/K 
for B || ab and B || c, respectively, corresponding to a 
— 2.8 and 0.59. The coherence lengths are calculated to 
be £ Q b = 7.2 nm and £ c = 1.5 nm. The Maki parame- 
ter a for B || ab is so large that theories developed for 
single-band superconductors would predict a first-order 
superconducting transition or the FFLO state. Further, 
the clear flattening at low temperatures of the B C 2(T) 
curve for B \\ ab indicates a strong paramagnetic limit- 
ing, and B C 2(0) — 12.7 T is close to a simple estimate 
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FIG. 2. (color online), (a) Superconducting resistive transi- 
tions in sample H2 at selected temperatures. The curves are 
vertically shifted for clarity. Hysteresis is visible at T — 0.03 
and 0.80 K but not at 1.1 K. The bottom curves are the field 
derivative dR/dB at 0.03 K. 
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FIG. 3. (color online) (a) Temperature dependence of the 
upper critical field in sample HI and (b) enlargement of a 
low-temperature part for B || ab. Upward and downward tri- 
angles correspond to up- and down-field data, respectively. 
The dashed curve for B \\ ab is a WHH fit in the whole tem- 
perature range, while the solid one is a fit in a temperature 
range T > 2 K. The solid curve for B \\ c is a WHH orbital 
limit curve. 
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FIG. 4. (color online) (a) Angle dependence of the upper crit- 
ical field in samples HI and H2. Upward and downward trian- 
gles correspond to up- and down-field data, respectively. The 
dashed and solid curves are calculated with the anisotropic 
Ginzburg-Landau theory. The former do not include the para- 
magnetic effect, while the latter do. 



of the paramagnetic critical field from T c , B po = 15.1 
T. A WWH fit22i to the a6-plane data in the whole tem- 
perature range gives the dashed line, which deviates up- 
wards around T = 2 K and then downwards below 1 
K [Fig. 3(b)]. If we use only data points above 2 K, 
we obtain an excellent fit down to 2 K as shown by the 
solid line. In either fit, B C 2(T) decreases with decreasing 
T below T -2 K. Within the WHH theory^ this sug- 
gests that the superconducting transition becomes first 
order in the low temperature region. Experimentally, 
the hysteresis is observed below ~1 K [Figs. 1(a) and 
2]. The phase boundary below ~1 K is nearly parallel 
to the T axis, and hence it follows from the Clapeyron 
equation dB/dT = -AS /AM that the entropy differ- 
ence between the superconducting and normal phases at 
the phase boundary is nearly zero. It would therefore be 
very difficult to see this phase boundary via heat capacity 
measurements 

We also note that the experimental i? C 2(0) of 8.0 T 
for B || c is larger than the orbital critical field B* 2 (0) 
— 6.4 T estimated from the initial slope. The solid line 
drawn for the c-axis data shows a B C 2(T) curve for the 
WHH orbital limit without the paramagnetic effect. The 
upward deviation of the experimental data indicates the 
importance of multi-band effects^ 

Figure 4 shows the angular dependences of B C 2 for 
samples HI and H2. Within the anisotropic Ginzburg- 
Landau (GL) theory B* 2 (9) = B* c2 {6 = O)/5(0) 

with 5(9) = \/cos 2 9 + e 2 sin 2 (9, where e = B* 2 {9 = 
0)/B* 2 (9 — 90°).— Neglecting the paramagnetic effect, 
namely identifying B* 2 with B c2 , these formulas are often 
used to describe the angular dependence of B c2 . How- 
ever, they fail in the present cases. The dashed curves 
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FIG. 5. (color online) (a) Magnetic torque in sample H2 mea- 
sured at 9 — 87°. The R(B) curves at the same angle are 
shown for comparison, (b) Magnetic torque for B \\ ab within 
experimental accuracy (~0.5°). The solid and dashed curves 
correspond to T = 0.03 and 1.5 K, respectively. Note that 
the vertical scale is roughly two orders-of-magnitudes smaller 
than that in (a). The R(B) curves at 0.03 K for B || ab are 
shown for comparison. Kinks appear in the torque curves at 
0.03 K as indicated by the grey vertical bars attached to the 
field-down curve. The kink structures may also be explicable 
as due to two-step jumps in the torque as illustrated with the 
field-up curves. The hysteresis in the torque up to the highest 
shown field may be ascribed to surface superconductivity (see 
text). 



are drawn using the experimental values of B* 2 (9 = 0) 
and e. They clearly deviate from the experimental data. 
When the paramagnetic effect is important, the angu- 
lar dependence of B c2 is modified: B c2 {9) = [B* 2 (9 = 
0) — aB^ 2 (9)]/S(9), where a is a parameter describing the 
strength of the paramagnetic effect X This gives excellent 
fits to the experimental data as shown by the solid curves, 
again confirming the presence of the strong paramagnetic 
effect. The fitted parameters are [B* 2 (9 = 0), a, e] — 
[11.6(4), 0.058(5), 0.17(4)] and [8.5(1), 0.067(2), 0.07(2)] 
for samples HI and H2, respectively, with B* 2 (9 — 0) and 
a -1 in Tesla. 

Finally, we show magnetic torque data measured on 
a small piece cut from sample H2. Figure 5(a) shows 



the torque measured for the field direction 9 = 87° over 
a wide field range. The U-shaped torque curves with 
a peak approximately at B c2 /2 can be explained in the 
framework of the anisotropic GL theory^ albeit the dif- 
ference between the field- up and -down curves, which is 
due to vortex pinning. The field-down curve indicates a 
weak peak effect in a field range ~9 to 10 T. The field-up 
and -down curves gradually merge above the peak effect, 
and the resistive transition occurs in this region as in- 
dicated by the R{B) curve measured at the same angle. 
The present torque curves bear a general resemblance to 
those reported for MgB2!^ 

Figure 5(b) shows the torque curves for B \\ ab within 
experimental accuracy (~ 0.5°) measured at T — 0.03 K 
(solid) and 1.52 K (dashed). The R{B) curve at T = 0.03 
K for B || ab is also shown for comparison. The torque 
curves at T = 0.03 K show kinks in the field region of the 
resistive transition as indicated by grey vertical bars for 
the field-down curve, and they approximately correspond 
to the onset of the resistive hysteresis, a kink in the R(B) 
curve, and the end of the hysteresis. The kink structures 
may also be interpreted as due to two successive jumps 
in the torque as illustrated for the field-up curves. The 
kink structures are not seen in the torque curves at T — 
1.52 K, where the resistive transition is not hysteretic. 
One might think that the observed jumps were too blunt 
for a first-order transition. However, we note that mag- 
netic torque measurements measure the component of the 
magnetization that is perpendicular to the field. Magne- 
tization jumps along the field direction could be much 
larger. Indeed, it has been reported that the afr-planc 
magnetization for B \\ ab in KFe2As2 exhibits an abrupt 
increase at B c2 at low temperatures^ 

The field-up and -down curves at T — 0.03 K do not 
merge in the shown field region, and the resistance curves 
are concave up to the highest shown field. These may be 
ascribed to surface superconductivity, which can persist 
up to B c3 ~ 1.7B C 2.— Although it is difficult to precisely 
determine B c3 in this sample because of drift in the mea- 
surement equipment (compare the two field-up curves 1 
and 3), the torque did not show clear irreversibility be- 
yond the drift when the field sweep direction was reversed 
at 17 T from up (curve 1) to down (curve 2). The surface 
superconductivity was previously reported for MgB2 ) 34 i 37 
for example. 

The magnetic torque for field directions close to B \\ 
c exhibits de Haas-van Alphen (dHvA) oscillations for 
fields above ~16 T, confirming the good quality of the 
sample. 38 Using the observed dHvA frequency of F w 
2000 T, we can estimate the carrier mean free path to be 
/ > 50 nm, which is much longer than £ a j, indicating that 
the sample is in the clean limit. 

Theoretically, the superconducting transition in simple 
paramagnetically limited superconductors becomes first 
order below 0.56T C £^ Since the first order region would 
shrink in realistic cases with the orbital effect, the present 
narrow hysteretic region (T < 0.1T C ) does not conflict 
with the theoretical prediction. Because the Maki param- 
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eter a for B || c is less than one, the paramagnetic effect 
becomes less important as the field is tilted from the ab 
plane, which explains the disappearance of the hysteresis 
for 8 < 80°. We also consider the FFLO states, since the 
samples are in the clean limit. An early theoretical study 
predicted a possible occurrence of the FFLO state for 
a > 1.84^ The maximum temperature of the FFLO re- 
gion was estimated to be 0.55T C for a = oo but was shown 
to decrease with decreasing a. With this model^ the 
transition from the mixed state to the FFLO state was 
found first order, while that from the FFLO to the nor- 
mal state second order. However, recent theories predict 
more complicated phase diagrams, in most cases with 
a first-order transition line in some parts of them, and 



whether the FFLO-to-normal transition is first or second 
order depends on parameters of models considered.— ~— 
Thus, the present observation of the resistive hysteresis 
does not contradict possibility of the FFLO state. 

IV. SUMMARY 

The temperature and field angle dependences of B C 2 
in Bao.o7Ko.93Fe2As2 indicate a strong paramagnetic ef- 
fect for B || ab, and the observed resistive hysteresis can 
be attributed to a first-order superconducting transition. 
The features observed in the magnetic torque data may 
be regarded as a sign of magnetization jumps. 
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clear peak at about F — 2000 T, which corresponds to the 
a frequency observed in KFe2As2 (The inset suggests 
the existence of another much lower frequency F ^200 T, 
but it is not resolved by Fourier transformation because 
of the too narrow field range.) The observation of the 
dHvA oscillations indicates that the Landau level spac- 
ing hu> c is comparable to the Landau level broadening F 
= U/2t due to impurity/defect scattering at these fields 
(B ^16 T). Using the observed dHvA frequency, we can 
estimate the carrier mean free path to be I = tvf 50 
nm, which is much longer than £ a (, indicating that the 
sample is in the clean limit. 



FIG. S 1. Fourier transform of dHvA torque oscillations at 6 
= 10° with background subtracted (inset). 



The magnetic torque for field directions close to B || c 
exhibits de Haas-van Alphen (dHvA) oscillations. The 
inset of Fig. 1 shows the oscillatory part of the torque 
for 8 = 10°. The Fourier transform (main panel) shows a 
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